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Abstract 

We present an integrability criterion for rational mappings based on two requirements. First, 
that a given point should have a unique preimage under the mapping and, second, that the 
spontaneously appearing singularities be confined to a few iteration steps. We present sev- 
eral examples of known integrable mappings that meet these requirements and, also, use our 
algorithm in order to derive new examples of integrable mappings. 



1. Introduction 

Integrability detectors are rare even for continuous systems. Recent progress in the integrability 
of discrete systems has spurred the activity in this direction, leading to interesting results of great 
variety. In order to characterize discrete integrability, Arnold has introduced and investigated 
the concept of complexity for mappings in a plane [1] . Arnold defines the complexity from the 
number of intersection points of a fixed curve with the image of a second given curve under 
the k th iteration of the mapping. For a polynomial mapping the growth of the number of 
intersection points is in general exponential in k. However for integrable mappings the growth 
is only polynomial in k. This result is included in a more general analysis presented by Veselov 
[2] discussing the dynamics of multiple-valued mappings (correspondences) and the growth of 
the number of different images (and preimages). Veselov, too, has linked integrability to "slow 
growth" . 

Our approach of discrete integrability is different from the above since it was based, 
essentially, on rational (rather than polynomial) mappings. For rational mappings, an important 
question is what happens whenever accidentally (i.e. depending on the initial conditions) a 
denominator vanishes, leading to a divergent mapping variable. In general one expects this 
singularity to propagate indefinitely under the mapping iterations, but it turns out that for 
integrable mappings these singularities disappear after a few steps. This observation has led to 
the proposal of the singularity confinement [3] criterion as detector of discrete integrability. Its 
efficiency has already been proven through the derivation of new integrable systems leading to 
the discovery of discrete Painleve equations [4]. 

It appears that the two notions of slow growth and confined singularities play an impor- 
tant role in the characterization of the integrable discrete systems. In what follows, we will try 
to present our approach which is based on both notions. We will, first, introduce the notion of 
preimage nonproliferation as well as the algorithm for its assessment. Based on the slow-growth 
principle, we claim that the number of preimages of a given point should not grow exponentially 
fast, which, when we consider mappings rather than general correspondences, can only mean a 
single preimage. Therefore, for the practical implementation we will present of this criterion, we 
will require that the inverse of the mapping be uniquely defined. The singularity confinement 
conjecture will also be extended in the following sections, essentially through the extension of 
the notion of singularity of a mapping. Apart from the singularity related to a divergence, we 
will consider as appearance of a singularity all the instances where the mapping accidentally 
loses some degrees of freedom. (The precise mechanism will become clear in Section 3). Con- 
finement of this singularity consists in the recovery of these lost degrees of freedom usually 
through the appearance of an indeterminate form like 0/0. While preimage nonproliferation is 
only a necessary condition for integrability, we conjecture that its combination with singularity 
confinement leads to a sufficient condition for integrability of discrete systems. 

In what follows, we will limit ourselves to rational explicit mappings, i.e. 

x'i = fi(x 1 ,x 2 ,...x N ) i = l,...N (1) 
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with rational /j's. Integrability in our sense means one of the following things: 

a) existence of a sufficient number of rational $fe(a;i, . . . Xn) = Cfc the values of which are 
invariant under the action of the mapping. 

b) linearizability of the mapping through a Cole-Hopf type transformation Xi = Pi/ Qi where- 
upon the mapping reduces to a linear one for the Pi's, Qi's. 

c) linearizability through a Lax pair. In this case, the mapping is the compatibility condition 
of a linear system of differential-difference, q-difference or pure difference equations. 

The above are not definitions but rather illustrations of the various types of integrability. 
It may well occur, as in the case of Quispel's mappings [5], that the existence of one invariant 
reduces the mapping to a correspondence of the form F(x, x') = that can be parametrized 
in terms of elliptic functions. In other cases, integration using the rational invariants may 
lead to some transcendental equation like the discrete Painleve ones. All of the above types 
of integrability have been encountered in the discrete systems that we have studied [6,7]. The 
reason for the above classification is to emphasize the parallel existing between the continuous 
and discrete cases. In the next section , we examine specific examples of mappings and formulate 
along the way our conjecture on preimage nonproliferation. 

2. Examples of integrable mappings and the preimage nonproliferation 
criterion 

Let us start with a very simple example of rational mapping, in which the growth of the number 
of preimages must be invoked. In [6], we studied the one-component, two-points mapping of 
the form: 

x' = f(x) (2) 
where / is rational. Singularity confinement considerations lead to 

f( x ) = a + Y — L__ ( 3 ) 

Y (x - v ; 

with integer Vk, provided that for all k, (3k ^ ce. Indeed, if x = /3& at some step, then x' 
diverges, x" = a and x'" is finite. So the mapping propagates without any further difficulty. 
However, if we consider the "backward" evolution, then (2) solved for x in terms of x' leads 
to multideterminacy and the number of preimages grows exponentially with the number of 
"backward" iterations. Indeed, the only mapping of the form (2-3) with no growth is just the 
homographic: 

, ax + b 

X' = : (4) 

cx + d v ' 

which is the discrete form of the Riccati equation. Thus, in this case, the argument of slow 

growth of the number of preimages of x is essential in deriving the form of the discrete Riccati 

equation. 

Another classical example in the domain of integrable mappings is the Quispel familly. In 
[5], Quispel and collaborators have shown that the mappings 

h(y) ~ f2(y)x 



x 



hiy) - h{y)x 



(5) 

_ 9i(x)-g 2 (x)y 

y = — — — — 

92{x) -9z{x)y 

are integrable, provided the /j's, g^s are specific quartic polynomials involving 18 parameters. 
We remark here that the mapping is "staggered", i.e. while x is defined in terms of (x,y), y 
is defined in terms of (x,y). It is precisely this staggered structure that allows one to define 
a unique preimage to (x,y). As Quispel has shown in [8], the mapping (5) is reversible which 
means that it can be written as a product of two involutions. It is not clear whether reversibility 
is a prerequisiste for integrability but though reversibility ensures that the preimage is unique, 
still there exist reversible systems that are not integrable. 

Reversible integrable mappings have also been considered by the Paris group [9] in their 
works based on the study of lattice spin and vertex models. They have shown that the trans- 
formations involved are in fact symmetries of the Yang-Baxter equations. These symmetries 
are constructed as the product of a pair of noncommuting involutions: thus the mapping is 
reversible and generically of infinite order. Still, it is interesting to investigate the mechanism 
for the nonproliferation of preimages in this case and present the algorithm that one should use. 
Let us illustratre this in the case of the mapping: 

, x + y - 2xy 2 
x — 



y(y - x) 

(6) 



, x + y — 2yx 2 
V x(x - y) 

The first step consists in considering the system of N equations x\ — fi(xt) = and eliminate 
successively all the x'^s but one. In the case of the mapping (6) the resultant in x after elimi- 
nating y (and vice-versa) is a fifth-degree polynomial in x (resp. y) with coefficients depending 
on x' and y' . Next we factorize this resultant. Preimage nonproliferation requires that only one 
factor depend on x',y', the other factors being associated to indeterminate forms 0/0. We find 
in the particular example the factors x 2 [x 2 — 1) (resp. y 2 {y 2 — 1)) and one last factor leading to 

x'-y' 



y ,2 + X l y l _ 2 

(7) 



y' 

y = - 



.12 _|_ X ,yl _ 2 

This is the typical situation for integrable rational mappings. The factorization of the resultant 
gives the unique inverse of the mapping along with particular values (here x = y = 0, ±1) 
corresponding to the indeterminate forms of the mappings. 

Veselov has studied the integrability of polynomial mappings and has shown [2] that the 
mapping x' = P(x,y),y' = Q(x,y) is integrable (in the sense that it possesses a nonconstant 
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polynomial integral $(x,y)) if there exists a polynomial change of coordinate variables trans- 
forming the mapping to triangular form: 



x 1 = ax + P{y) 

V' = PV + 1 (8) 

for polynomial P. Moreover he has shown that in this case the complexity of the mapping 
is bounded. The important feature in (8) is the fact that the equation for y' is linear. Thus 
the inversion of (8) is straightforward. Thanks to the triangular form the integration of (8) is 
reduced to the solution of two affine mappings, first for y and then for x. 

One more interesting illustration of the preimage nonproliferation algorithm is provided 
by the discrete Painleve equations that we derived in [10] and which are not of Quispel form. 
We have found there that the mappings: 

/ = xy(a(y+ 1) - xy 2 ) 
a(y + l? 

(9) 

/ = a{y + l){xy 2 -(y + l)(a- zy) 
V ~ (a(y+l)-xy 2 ) 2 

where a = cnst. and z linear in the lattice variable, is a discrete form of the Pi equation. In 
order to check the preimage nonproliferation, we eliminate x (or y) from (9) and factorize the 
resultant. We find as expected factors related to pathological points x = 0, y = 0, y = — 1, and 
a unique inverse that reads: 



x = 



x'(x'y' + z)[(x'y' + z) 2 + a(y' + l)(z - x')\ 
a(y' + l) 2 

(10) 

= a{z-x>){yi + 1) 
(x'y' + z) 2 

Thus in this example, too, as in all previous ones, integrability is related to non-growth of the 
number of preimages. 

2. Extending the singularity confinement criterion 

In the previous section, we encountered several examples of integrable mappings, all of which 
satisfied the "no-growth" property. Here, we will apply the preimage nonproliferation criterion 
in order to construct explicitely integrable mappings. However since this criterion furnishes only 
a necessary condition for integrability, we will supplement it by singularity confinement, our 
conjecture being that the combination of the two criteria is sufficient for discrete integrability. 

Applying the preimage nonproliferation algorithm to a general mapping can easily lead to 
untractable calculations. If, however, there are not too many free parameters in the mapping 
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the implementation of the criterion is straightforward. In what follows, we will limit ourselves 
to simple two-component, two-point mappings of the form 

/ Qi(x,y) , Qi(x,y) 

Q(x,y) V Q(x,y) [ ] 

where Q, Q\ and Q 2 are quadratic polynomials in x and y. By applying a general linear 
transformation on this mapping we can reduce the (common) denominator to one of the two 
canonical forms Q = xy — 1 or x 2 — y (or any of the degenerate forms Q = xy, x 2 — 1 or x 2 ). 
Let us start with the mapping h: 

, _ a2ox 2 + anxy + ao 2 y 2 + aipx + a 01 y + a o 
xy — 1 

, b 20 x 2 + buxy + b 02 y 2 + b w x + b 01 y + b 00 

V = ; (12.6) 

xy-1 

One assumption that we will introduce here is 002 = since it leads to a great simplification 
of the calculations. As a first step in the preimage nonproliferation algorithm we eliminate y 
between (12.a-b) for given x' and y' and obtain a resultant that is a fourth degree polynomial 
in x. We demand that three of the roots be independent of x' , y' and denote them by x±, x 2 
and X3. We then demand that whenever x = Xi, y = i = 1,2,3 both x' and y' have the 

indeterminate form 0/0. Calling: 

S = x\ + x 2 + x 3 
P = x t x 2 + x 2 x 3 + X3X1 (13) 
II = X1X2X3 
we find condition for a unique preimage: 

a i = -a 2 oII 

aoo = 020P - an (14) 
610 = -620S - 6 02 /n 

b 01 = -b 20 u - 6 02 p/n 
&oo = 6 20 p - fen + &02X/11 

With (14) the mapping satisfies the preimage nonproliferation requirement. This is not, however, 
sufficient for integrability. What we must also demand is that the mapping have confined singu- 
larities. The simplest kind of singularity is whenever the denominator vanishes, i.e. Q(x, y) = 
or in our example y = 1/x. However, since in the present case the numerators Qi, Q 2 are 
also quadratic, the singularity is confined in one step: Q(x,y) = leads to diverging x', y' and 
because the degrees of numerators and (common) denominator are equal this leads to finite x" 
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and y" . So the study of this singularity does not introduce any constraint on the mapping. 
But the vanishing of the denominators is not the only singularity of the mapping: a subtler 
singularity may exist. 

Normally for a general ./V-component mapping, N free parameters, introduced by the 
initial conditions, must be present at every step. Now, it may happen that at some iteration 
one (or more) degress of freedom be lost. The condition for this to occur is that the Jacobian 
of (x[, x' 2 , ■ ■ ■ x' N ) with respect to (xi,X2, ■ ■ ■ xn) vanishes. For a general mapping x- = fi(xk) 
this reads: 



J 



dx j 


dx[ 


dx' ± 


dxi 


dx 2 


dx jv 


9x' 2 


9x' 2 


9x' 2 


dxi 


dx 2 


dx jv 


dx' N 


dx' N 


dx' N 


dxi 


dx 2 


dx jv 



(15) 



How can this singularity be confined? By this we mean that the mapping must recover the lost 
degree of freedom. For rational mapping of the kind we are considering, this can be realized 
if some of the mapping's variables assume an indeterminate form 0/0. In that case new free 
parameters can be introduced and the mapping recovers its full dimensionality. 

Let us apply this criterion to the mapping h. The Jacobian readily factorizes and we 
obtain three factors: 

x + Xi - E + yU/xi = i = 1,2,3 (16) 

Thus whenever (16) is satisfied a singularity appears (in the sense of the loss of one degree of 
freedom). For the confinement of this singularity (at the x", y" level) we must have x'y' — 1 = 
and the numerators of both x" and y" must vanish. First we supplement the condition for the 
vanishing denominator x'y 1 — 1. This leads to a number of equations that, in fact, specify fully 
the remaining a, b coefficients and moreover put a constraint on x\,x 2 ,xs- The latter can be 
written (up to an odd permutation of xx,x 2 ,xs) as: 



3xix 2 x 3 



2 i 2 i 2 
X^3 + X 2 X\ + X 3 X 2 



(17) 



For the a, b we obtain: 



O20 = - 
b 2 



-n 
i 

an 
an 



(18) 



ho = n 

&n = 77 

II an 

In fact the simplest way to parametrize equations (14,17,18) is to take an = a as basic parame- 
ter. Introducing one further parameter uj we can express Xi,x 2 ,x% as: X\ = —auj, x 2 = a(l+l/u) 
and X3 = a/(l + uj). It turns out that once conditions (14,17,18) are implemented the numera- 
tors of both x" and y" automatically vanish. Thus the mapping h is singularity confining and 
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according to our conjecture it should be integrable. This is indeed the case and one invariant 
can easily be found. It reads: 

fl X + Xi x-x VlL/Xi ={ - irK (19) 
i=i 1 

i.e. the product on the lhs, instead of being strictly constant, alternates sign between even 
and odd iterations. One should, in principle, take the square of the lhs in order to find a true 
constant. A closer inspection of the mapping (motivated by the form of the invariant) reveals 
an even simpler structure: the mapping is periodic with period six, i.e. h 6 = I. 

Finaly, the mapping can be cast in a much simpler form if one use the scaling freedom in 
order to reduce the number of the parameters from two to one. We shall not enter into these 
details but we just give the final result: 

, —x 2 + xy + ax — y + 2 — a 

x' = (20.o) 

xy — 1 

, = -x 2 + (2 - a)xy + y 2 + (1 + a)x + (a - A)y + 1 - a 
V xy-1 

In an analogous way we can treat the mapping p: 

, a 20 x 2 + a n xy + a 02 y 2 + a 10 x + a 01 y + a 00 , 01 v 

x = 5 (21. a) 

x z — y 

, b 20 x 2 + buxy + b 02 y 2 + b w x + b 01 y + fe o , 01 
y = 21.6) 

x l — y 

with 002 = 0. As in the previous case, we can ask that the resultant of the elimination of y 
between (21.a-b) for given x' , y' (which is quartic in x), possess three roots independent on x' , 
y', namely X\, x 2 and X3, corresponding to 0/0 indeterminacies. In a second step, the singularity 
confinement can be implemented in a perfect parallel to the case of the mapping h, leading to: 

*' = x2S " X r XF + n (22.a) 
x A — y 

, = x 2 E 2 - 2xj/S + y 2 - (SP + n)x + yP + Sn 

y x 2 — y 

with S, P, II given by (13). As in the previous case, the freedom of transformations can be used 
in order to simplify this mapping. Finally only one parameter remains and the mapping reads: 

x' = x ^-y-^ (23 .a) 
x z — y 

y> = ( x -vX*-y-P) (23.6) 
x 2 — y 

This mapping is indeed integrable, but in a trivial way: it is just an involution, p 2 = I. Still, 
the important point here is that the conjecture concerning the integrability of mappings that 
have nonproliferating preimages and confined singularities is once more satisfied. 
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Conclusion 

In the preceding sections, we have presented in detail the preimage nonproliferation criterion, 
which, we conjecture, is a necessary condition for the integrability of rational mappings. Based 
on the "slow-growth" principle, this criterion consists in requiring that, in order to be a candidate 
for integrability, a rational mapping possess a unique inverse. Thus the number of preimages 
of a given point through the mapping does not grow with the number of iterations (while 
an exponential increase is, generically, expected). Since this criterion offers only necessary 
conditions it cannot predict integrability but can be used as a fast "screening" procedure. The 
successful candidates can then be tested for singularity confinement, which is more stringent but 
of more difficult implementation. The combination of the two criteria (preimage nonproliferation 
and singularity confinement) we conjecture to be an integrability predictor of the same efficiency 
as the Painleve method for continuous systems. Several results exist already based on this 
approach and we expect the extension of the singularity confinement presented here to further 
widen its range of applications. 
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